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Abstract 

The neutron-drip-line nucleus 23 N is investigated in a three-body model consisting of a 21 N core 
and two valence neutrons. Using the Faddeev formalism with the realistic neutron-neutron potential 
and the neutron-core potential, we calculate the ground-state properties of 23 N including its two- 
neutron separation energy, and obtain good agreement with the experiments. We also find an 
excited 23 N state with a shallow two-neutron separation energy at about 0.18 MeV. By evaluating 
the root-mean-square matter radii, the average distances between the two valence neutrons, and 
the average distances from the core to the center-of-mass of the valence-neutron pair, we show that 
the excited state of 23 N has a distinct halo structure. Through calculating the correlation density 
distributions of the 23 N three-body system in configuration space, we find that the excited state 
of 23 N has a triangular shape, which is similar to but much more extended than the ground state. 
This scaling symmetry between the ground and excited states indicate the existence of an Efimov 
state in 23 N. 
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I. INTRODUCTION 


The properties of neutron-rich nuclei close to the drip-line are crucial for understanding 
fundamental mechanism governing the stability and formation of nuclei. In recent years, 
tremendous experimental and theoretical efforts have been devoted to the investigation of 
neutron-rich nuclei [1-5| . The discoveries of light nuclei with halo structures have attracted 
great attention because of their exotic properties 6-tlQ]. and their important role in the big 
bang nucleosysnthesis jll]. For the nitrogen isotopes, many studies have focused on neutron- 


rich isotopes with mass number ranging from A = 17 to A = 22 


12 


17]. However, 23 N has 


rarely been investigated since its first production in 1998 jl8|, [19|. As a newly synthesized 
drip-line nucleus, many physical properties of 23 N such as the energy spectrum and the 
structure configuration have not been observed or predicted yet except for the one-neutron 
separation energy S n . S n s in 22 N and 23 N are respectively 1.28 MeV and 1.79 MeV, which are 
much smaller than that in 21 N of 4.59 MeV [20]. Therefore, in a first-order approximation, 
one can treat the 23 N nucleus as a three-body system composed by an inert 21 N core, which 
is tightly bound, and two valance neutrons moving at a relatively large distance away from 
the core. 

Various few-body methods have been adopted to calculate different three-body quantum 
systems, such as the Faddeev formalism |2lJ, the equivalent two-body method 22|, |23J, the 


Green’s Func 
function [25 


ion Monte Carlo 


24], and the THSR (Tohsaki-Horiuchi-Schuck-Ropke) wave 


26]. Among these different approaches, the Faddeev formalism has a specific 


advantage that it can also describe the three-body mechanism in a heavier system, which is 


normally encoded in microscopic calculations, in a computationally simpler way 2l|. Various 
neutron-rich nuclei, including 6,8 He, n Li, 12,li Be, 1 'B, and 22 C have been investigated with 


the Faddeev formalism 


27 


This approach has also been applied recently to proton- 


rich nuclei, such as 1 'Ne, 18 Ne, and 28 S [37|, [38]. In this work, we have applied the Faddeev 
equations to the 21 N+n+n three-body system. By solving these equations, we obtain various 
properties of the 23 N nucleus, which are essential for understanding the halo structure inside 
the nucleus. 

In Sec. HU we provide a brief introduction to the Faddeev equations and the two-body 
interactions adopted in our calculations. In Sec. IIII1 we present the numerical results of the 
three-body calculation and analyze the halo structure of the ground state and our newly 
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discovered excited state in 23 N. At last, the conclusion is given in Sec. IIVI 


II. FORMALISM 

In the Faddeev formalism, the full wave function T JiU for a three-body system can be 
decomposed into three components with respect to different sets of Jacobi coordinates 


V JM = *i M {x 1 , 2 / 1 ) + ^ J 2 M (x2,y 2 ) + (+ 3 - //:’.)• 




JM( 


,JM / 


( 1 ) 


where ^ J i AI is the z-th component of the full wave function with total angular momentum 


J and its z-component M. The i-th set o 
component is shown in Fig. Q] for i = 1,2,3 


re Jacobi coordinates corresponding to this 
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FIG. 1. Three sets of Jacobi coordinates in the Faddeev formalism for the core+n+n three-body 
system. 

The coupled-channel Faddeev equations for the core+n+n system are expressed as 


V\ - E)m{ M = + T 3 


,jm __ 

,JM 


(Ti + V 1 - E)^{ M = —Vi(^> j 2 m 

(T 2 + v 2 - E)^i M = 

(t 3 +1/ 3 - e)^>i m = -v^{ M + t J 2 M ) 

where Ti = T X i + T y i is the relative kinetic energy in z-th coordinate set. V* = Vjk represents 
the effective two-body nuclear force between particles j and k. The indices (z, j, k ) = (1, 2, 3) 


3 ) 
JM) 

JM\ 


( 2 ) 


run through cyclic permutation 


2ii], 
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We then introduce the Jacobi coordinates, which include the hyper radius p and the 


hyper-angle 6U They are related to the Jacobi coordinates by 


21 


2 2,2 
P =Xi+Vi 


,Xi 


Qi = arctan( —). 

Vi 


(3) 


One should notice that the hyper-radius is the same in different sets of Jacobi representation, 
but the hyper-angles are different. To solve the Faddeev equations, we adopt the hyperspher- 
ical harmonic method. Using this approach, one can conveniently separate the hyper-angle 
and hyper-radial dependencies of the wave function. The hyperspherical decomposition of 
the i-th component of the wave function is defined as 29|, l37| 


= p-** E E x&fcWtfrw i*=“<>, 


(4) 


where cq = { l X i , l y i, Li, Sj, Sk, S X i}Ji and |i : ati) represent the quantum numbers in the spe¬ 
cific partial-wave channel occupied by the three-body system. The hyper-angular function 
(j)^ lyi (Oi), which are the eigen-solution of the hyper-angular equation, is expanded in terms 
of Jacobi polynomials as 


<j>^: yi (Qi) = N l £ yi (sin Qi) L ™ (cosQi) Lyi P l ™+ 1/2 ’ lyi+1/2 (cos2 Q { ). 


(5) 


Here Pn": z+1 ^ 2 ' lyz+1 ^ 2 (cos20i) denotes the Jacobi polynomial with n* = 0,1,2,---. imp:”' is 
the normalization coefficient, and Ki = l XA + l yt + 2n t indicates the hyper-angular momentum 
with respect to the corresponding Jacobi polynomial. 

After substituting Eq. (J5|) into Eq. (QJ, the hyper-radial function Xq^k, (p) an d the hyper- 
angular function (p '^" y% (Qi ) are obtained using the numerical program to solve the Faddeev 


i\T^xi^yi 
ly Ki 


equations 


2111 - 


In the calculations of 23 N as the 21 N+n+n system, the core-neutron and the neutron- 
neutron interactions can be determined phenomenologically by fixing the experimental data. 
In this paper, we adopt the well-known GPT potential neglecting only the spin-spin term 
for the n-n interaction and keeping the central, tensor and spin-orbit parts. This potential 
provides good fits to the low-energy properties of the low-energy n-n scattering (tL, 40]. For 


the core-neutron interaction, we adopt a Woods-Saxon form 


m 


U 




U n _ c 


Un 


r = 


+ 


14 


exp( 


r—r o v 
a 


l + exp(^) ra (1 + exp( r ^ m )) : 


;L'S, 


(6) 
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where r 0 = 1.25A 1 / 3 fm and a = 0.65 fm 32, [4l|. Vo is the depth of the Woods-Saxon 
potential, and I4 0 represents the strength of the spin-orbit coupling. Vo and V so can be 
determined by fixing the binding energies of the core-neutron two-body subsystem. We 
then apply the super symmetric transformation to this interaction to eliminate the spurious 
bound states which are forbidden by the Pauli principle. 


III. NUMERICAL RESULTS 


We discuss in this section the numerical results for calculating 23 N in a three-body model. 
We firstly discuss the physical model for the core and explain our choices of parameterization 
for the effective n- 21 N interaction, which are suitable to reproduce the experimental data. 
Then we solve the Faddeev equations and calculate the neutron-separation energies and 
configuration-space wave functions of 23 N. These observables can be utilized to analyze the 
halo structure of 23 N. 

To construct the effective n- 21 N interaction, information about the shell-model occupa¬ 
tions of the 21 N core and valance neutrons, which can be determined from experiments, 
needs to be known. Here we make two assumptions in the shell-model description. Firstly, 
we assume that the neutrons inside the 21 N core occupy the (Osi/2) 2 , (OJO3/2) 4 , (Opi/2) 2 and 
(0d 5 / 2 ) 6 orbits, which are forbidden to be occupied by the valence neutrons due to the Pauli 
principle |3l|. With the lack of the experimental information about the excited states in 


21 


N, we simply neglect these excited states in our calculations as a first-order approximation. 


Similar assumptions are also made in the work by other groups, e.g., in Refs 


• HQ, 


42]- 


With these limitations, we assume that the valence neutron occupies the lsi /2 state in the 
n-core sub-system of 22 N. 

Using the above assumptions and the Wood-Saxon potential form, i.e., Eq. (] 6 ]) , we de¬ 
termine the core-n interaction by reproducing two experimental results. The first is the 
one-neutron separation energy of 22 N, corresponding to the binding of the lsi /2 orbit, i.e., 
S n = e(lsi/ 2 ) = 1.28 ± 0.21 MeV. We also take into account the one-neutron separation 
energy of 21 N, which is related to the binding of the 0d 5 / 2 orbit, i.e., e(0d 5 / 2 ) = 4.59 ± 0.11 
MeV. The experimental error of e(lsi/ 2 ) is considered in our calculations; while the error 
bar of e(0d 5 / 2 ) is neglected due to its relatively small size. By reproducing the upper bound, 
lower bound and mean value of e(lsi/ 2 ), and also fixing the mean value of €( 0 ^ 5 / 2 ), we de- 
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termine three sets of parameters of the core-n interactions and list them in Table [U The 
corresponding root-mean-square (r.m.s.) radii of the core, the valence neutron and total 
matter in the 21 N-n two-body system are calculated with respect to different parameteriza- 
tions and shown in Table [TT1 The neutron density distributions of 22 N are shown in Fig. [21 
where the probability density of the neutrons inside core and of the last neutron occupying 
the Is 1/2 state, are shown respectively. We find that the last neutron has a much more 
extended matter density distribution than the neutrons inside the core do. This clearly 
indicates a halo structure in the 21 N-n system. 

TABLE I. Parameters of the effective 21 N-n interaction, e is the single-particle energy of the 21 N-n 
system. The upper bound (1.49 MeV), the mean value (1.28 MeV), and the lower bound (1-07) 
MeV of S n in 22 N are adopted for sets A, B, and C, respectively. 




^0 

a 

F 0 

F so 

e (0 s l/2) 

e(0P3/2) 

e(0pi/ 2 ) 

e(Odj5/2 ) 

e (l s l/2) 



(fm) 

(fm) 

(MeV) (MeV.fm 2 ) 

(MeV) 

(MeV) 

(MeV) 

(MeV) 

(MeV) 

Set 

A 

3.45 

0.60 

-41.45 

-46.45 

-25.382 

-14.939 

-9.026 

-4.590 

-1.490 

Set 

B 

3.45 

0.60 

-40.77 

-50.51 

-24.813 

-14.647 

-8.224 

-4.590 

-1.280 

Set 

C 

3.45 

0.60 

-40.02 

-55.00 

-24.170 

-14.316 

-7.333 

-4.590 

-1.070 


TABLE II. The calculated r.m.s. radii of the core R^, the valence neutron /?*( 1 and total matter 
in 22 N. S n denotes the corresponding one-neutron separation energy. 




Sn 

nth 

rC c 

■nth 

n 

nth 

JX m 



(MeV) 

(fm) 

(fm) 

(fm) 

Set 

A 

1.490 

3.21 

5.25 

3.53 

Set 

B 

1.280 

3.22 

5.47 

3.58 

Set 

C 

1.070 

3.24 

5.75 

3.65 


With the appropriate parameterization for the effective 21 N-n interaction obtained above, 
we calculate the bound-state observables of 23 N in the 21 N+n+n three-body model using the 
Faddeev formalism. In these calculations, we use the cutoff K max =20 in the hyperangular- 
momentum expansions in Eq [4] to provide a proper numerical accuracy. Comparing with 
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FIG. 2. Density distributions of neutrons in the subsystem 22 N. The solid curve is the density 
distribution of the core neutrons, and the dashed curve is the density distribution of the last 
neutron. 


the results of K max =lQ, the difference between the two-body separation energies is about 
4%, which is much smaller than the relative error of about 30% due to different choices 
of parameters. Thus, this cutoff K max —20 is adequate for our Faddeev calculations of the 
three-body system. We calculate the two-neutron separation energies S 2n and r.rn.s. matter 
radii of the ground and excited states of 23 N using different sets of parameterization, and list 
the results in Table UTT1 We observe that both the the ground- and excited-state S 2n decrease 
when the depth of the 21 N-n potential increases. We also find that the r.rn.s. matter radii 
of the two states increase respectively with the decrease of S 2n in the corresponding states. 
This correlation which show that 23 N is a standard halo nucleus, as discussed in 43]. Our 
obtained two-neutron separation energy of the \ ground state, S 2n = 3.6 ±0.5 MeV, agrees 
well with the experimental result, i.e., S 2n = 3.07±0.31 MeV 20]- The r.rn.s. matter radius 
of the ground-state 23 N is about 3.0 fm, which corresponds to the size of a stable nucleus 
with mass number A & 27. This radius is much smaller than 22 C, whose r.rn.s. radius is 
about 5.4 ±0.9 fm 441 ]. Therefore, although loosely-bound compared with stable nuclei, the 


halo structure of the 23 N ground state may not be as clear as its isotone halo nucleus 22 C. 
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TABLE III. The two-neutron separation energies and r.m.s. matter radii r m of the ground 
and excited states of 23 N, calculated with three sets of parametrization and from the experiments. 
The superscript * denotes the excited state. 



S 2 n 

(MeV) 

r m 

(fm) 

Q* 

°2 n 

(MeV) 

1 m 

(fm) 

Set A 

4.13 

2.969 

0.315 

4.272 

Set B 

3.64 

2.985 

0.185 

4.358 

Set C 

3.13 

3.004 

0.069 

4.476 

Exp [20] 3.07T0.31 

- 

- 

- 


The d excited state of 23 N, which has not been discovered in experiments yet, is predicted 


in this work to have an extremely shallow two-neutron separation energy S% n . With different 
parameterization, we obtain S£ n to be in the region of 0.069-0.315 MeV, with a mean value 
of 0.185 MeV. The r.m.s. radii of this excited state is about 4.3 fm, which corresponds 
to the size of a stable nucleus with A ps 80. Therefore, the excited state of 23 N can be 
unambiguously described as a giant halo state. 

To illustrate the structure of 23 N, we calculate the average distances between the two 
valence neutrons r nn and the average distances from the core to the center-of-mass of the 
valence-neutron pair rc( nn ) in the 21 N+n+n system. The results in Table HV1 indicate that 
both ground and excited states of 23 N have triangular shapes. By calculating the ratios of 
r nn and rc(nn) respectively in the ground and excited states, we find that 


^C(nn ) 


C(nn) 


1.9 . 


(7) 


The same ratio obtained in both states suggests that these two states have similar geometric 
structures and are mainly different by a spatially discrete scaling factor. 

Furthermore, we also analyze the correlation density distributions of the ground and 
excited states of 23 N in configuration space. These quantities are evaluated in the Jacobi- 
coordinate representation with 21 N as the spectator particle, he., in the third diagram in 
Fig. m The spatial correlation density distribution is calculated as 33] 


P(r nn ,r c ,(nn)) = x z y 2 I \^- M (x,y)\' z dQ Xi dQ 


Vi- 


( 8 ) 










TABLE IV. The average distances between the two valence neutrons r nn and the average distances 
from the core to the center-of-mass of the valence-neutron pair rc< nn ) i n the 23 N ground and excited 
states The superscript * denotes the excited state. 




Tnn 

^C(nn) 

r* 

nn 

r C(nn) 



(fm) 

(fm) 

(fm) 

(fill) 

Set 

A 

6.415 

3.357 

14.600 

7.640 

Set 

B 

6.681 

3.496 

16.326 

8.543 

Set 

C 

6.878 

3.599 

17.059 

8.927 


The correlation densities for the ground and excited states are shown in Figs. [3] and [4], 
respectively. On the one hand, we find that the ground state has the largest probability 
density when the two neutrons are at a distance of about 3.5 fm from the core. On the 
other hand, the excited state has the largest probability density when the two neutrons 
are far away from the core with a distance of about 8.5 fm, and are separated from each 
other with a distance of about 16 fm. The spatial separation in the excited state of 23 N 
supports the halo structure suggested above. Furthermore, in Fig. [2 a second peak with 
a smaller amplitude, which mainly comes from the (Oaf) 2 component, also appears in the 
excited state. The occupancies of the valence neutrons are also calculated. For the ground 
state, the occupancies of the two valence neutrons in the (Is) 2 and (Oh) 2 states are 95% and 
5%, respectively. For the excited state, the occupancies of the two valence neutrons in the 
(Is) 2 and (Od) 2 states are 77% and 23%, respectively. 

Since the excited state has a much shallower two-neutron separation energy, a much larger 
r.m.s. radius, and much more extended spatial distribution than the ground state, one may 


be able to describe the excited state as an Efimov state in the halo nucleus 123 


45 


46]. 


In fact, one can utilize the geometrical simi 
configurations to distinguish Efimov states 


arity between the ground- and excited-state 


47]: One of the magic properties of Efimov 


physics is that two consecutive Efimov states can be related to each other by a discrete 

n 

spatial scaling factor (47J. In 23 N, the dominant parts of the spatial density distribution 
of the ground and excited states (Figs. [3] and [4j) indicate that the configurations of 23 N in 
these two states have very similar shapes. Moreover, the same ratio between r nn and rc(nn) 
in the ground and excited states (Eq. (|7])) suggests that a discrete scaling symmetry exists 
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FIG. 3. Contour diagram for the spatial correlation density distribution of the ground state of 23 N 
with parameters of Set B. (Color online) 



FIG. 4. Contour diagram for the spatial correlation density distribution of the \ excited state of 
23 N with parameters of Set B. (Color online) 


10 

















between the two states. These features are highly consistent with the explanation of Efimov 
states. However, further theoretical or experimental investigations are needed to confirm 
the excited state in 23 N. 

IV. CONCLUSION AND DISCUSSIONS 

We have investigated the properties of 23 N as a three-body system consisting of an inert 
21 N core and two valence neutrons. We apply the Woods-Saxon potential form to represent 
the effective 21 N-n interaction, and reproduce the basic characteristics of the 22 N. By solving 
the Faddeev equations of the 21 N+n+n three-body system, implemented with the 21 N-n and 
n-n potentials, we calculate the bound-state observables in the 23 N ground state, where we 
obtain the two-neutron separation energy S- 2 n consistent with the experimental result. We 
also discover an excited 23 N with a much shallower two-neutron separation energy of 0.18 
MeV. We also calculate the r.m.s. matter radii for the two states. The obtained results 
suggest a relatively small halo structure for the ground state. On the other hand, a more 
extended distribution of the valance neutrons is unveiled in the excited state, which indicates 
an distinct halo structure. The average distances between the valence neutrons, the average 
distances from the core to the center-of-mass of the valence-neutron pair, and the spatial 
correlation density distributions of the two states in 23 N are also calculated to investigate 
the geometric structure of the three-body system. These results indicate that the two states 
have very similar triangular shapes and can be related to each other by a discrete scaling 
symmetry. These features suggest that the excited state of 23 N can be an Efimov state. 
Our calculations of 23 N with the Faddeev formalism enable a new exploration toward the 
neutron drip line, and urge further experimental investigation of the 23 N ground state and 
the potential discovery of the 23 N excited state. 
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